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ABSTRACT 

The question of optimizing nozzle contours for micro-thrust 
rockets led to the design, construction, and testing of a low-density 
gas dynamic facility. The primary objective was to investigate the 
mass flow rates of a gas through various profiles in the slip and 
transition flow regimes at high pressure ratios. 

An initial test was conducted with an orifice as the test pro- 
file. The results showed that the facility can be used to investigate 
mass flow rates from the threshold of the free-molecule, through 
the transition and slip, to the continuum regimes. These results 
compare favorably with those of two previous investigators, and 
asymptotically approach the theoretical continuum and free-molecule 
limits. The ratio of mass flow rate to theoretical free-molecule 
mass flow rate is shown to transition smoothly from one theoretical 
Ties the other. A local maximum may occur inthis ratio in the 
slip regime,and the attainment of the theoretical free-molecule limit 


appears to occur more slowly than expected. 
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LIST OF SYMBOLS 

molecular BolGeion TAT etee: 

mean molecular speed 

aperture diameter 

transfer ratio of momentum 

gravitational constant 

height difference 

constant 

mass flow rate 

molecule number density 

pressure 

cylindrical or spherical radius 

axial velocity 

volume 

axial coordinate 

area 

discharge coefficient 

collision frequency per mole 

perimeter 

Boltzmann constant 

Knudsen number 

molecular weight 

Mach number 


number of molecules striking a unit area of wall per unit 
time 


pressure. 
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volumetric flow rate 

specific gas constant 

Reynolds number 

temperature 

volume 

characteristic velocity 

empirical coefficient 

Orifice discharge coefficient 

ratio of specific heats 

coefficient of slip 

spherical elevation angle or conical half angle 
mean free path | 
viscosity 

kinematic viscosity 

density 

spherical azimuth angle 

dimensionless flow function 


Clausing factor 


stagnation chamber or upstream 
exhaust chamber or dowmetream 
float 

iree-molecule 

free-molecule 


fluid 
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I. INTRODUCTION 
The flow of gases can be roughly divided into four regimes 
characterized by the Knudsen number (Kn), the ratio of molecular 
mean free path to a characteristic length in the flow field. These 


POUm ene ClImMes are: 


continuum flow ~ Kn <0. 01 
slip flow ~ 0.01 < Kn<0.1 
transition flow ~ 01 <Kn<i 
free-molecule flow ~ 5 <atar) 


The Navier-Stokes equations with appropriate boundary 
conditions describe the flows in the continuum and slip flow regimes 
while the kinetic theory of gases must be used in the transition and 
free-molecular regions. 

Ber flows through orifices, nozzles,and short tubes well 
into the continuum regime, the boundary layers on the walls usually 
represent only a very small portion of the flow field. These boundary 
layers become more and more dominant as the Knudsen number in- 
.creases, until they occupy the entire flow field at large Knudsen 
numbers. In the slip regime, the boundary layers no longer have 
zero velocity at the walls. In free-molecule flow,collisions between 
gas molecules can be neglected compared with collisions between gas 
particles and the body. In such cases, the fluxes of incident and re- 
flected particles can be treated separately. In the transition regime 
where the mean free path and the characteristic dimension of the 
body are oi the same order of magnitude, both intermolecular colli- 


sions and molecule-surface interactions must be considered. In 
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order to obtain experimental data for a better understanding of the 
transition from continuum to free-molecule flow, the present low- 
density gas dynamic facility was designed and built. 

The initial interest in this area of low-density gas flow 
was generated by the problem of nozzle contours for micro-thrust 
rocket motors. In these rockets with extremely low thrust, the ex- 
haust gases are in the flow regime between continuum and free- 
molecule. The optimization of these nozzle contours depends upon 
having some basic knowledge of gas flow in the slip and transition 
regimes. It is hoped that this low-density gas dynamic facility will 
supply some basic experimental data about siip and transition flow 


through various contours. 
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Ti APPARATUS 
General Description 

The low-density gas dynamic facility is a continuous flow, 
open-circuit type designed to operate with stagnation chamber pres- 
sures of from lu to $ atm and exhaust chamber pressures to are 
at least a 100:1 pressure ratio across the test section. 

Figure lis a photograph of the facility; subsystems are 
shown schematically in Figures 2-4. The pumping system consists 
of a 300-cfm rotary pump, Stokes Model 412-H, in parallel with a 
2-in. metal-oil fractionating diffusion pump, CVC Model PMCS-2C, 
backed by a 7-cfm rotary pump, Cenco Hypervac 25 (Figure 4). The 
pumping characteristics as supplied by the manufacturers are shown 
in Figures 5 and 6. 

The vacuum vessel (Figure 3) is a cylindrical steel ena 
ber of about 40-cm internal diameter and 90-cm length. The tank is 
divided into the stagnation and exhaust chambers by a central bulk- 
head which contains a cutout for mounting a machined test profile 
plate. For rapid pump-down, a bypass is incorporated around the 
test section. The tank has provisions for three stagnation-chamber 
pressure taps. The exhaust chamber pressure line enters the 
chamber through the 4-in. outlet and extends about half way into the 
chamber. The tip of the pressure line is off tO the side of, and 
turned at a right angle to, a direct line between the test section and 
-the outlet. Access to the test section is gained through an O-ring- 


sealed plate in the inlet side of the stagnation chamber. This plate 
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incorporates the inlet line and one pressure tap. 


Test Gas Feed and Measurement 

The test gas is fed from a standard gas bottle through 
first, a 25- to 650-psig high-pressure regulator, Matheson Model 2, 
and then, a 3- to 15-in. water-column low pressure regulator, 
Matheson Model 70 B, before passing to a liquid-nitrogen cold trap. 
From the cold trap, the gas flows through a constant-temperature 
water bath to the flow measurement system. At the inlet to the flow 
measurement system, there are both pressure and temperature taps. 
The pressure line is connected to a 0- to 15-psia diai gauge, Burton 
Model 2311-001 and a 0- to 800-mm mercury manometer, Wallace 
and Tiernan Model FA 135, and the temperature tap houses a pre- 
cision mercury thermometer, ASTM-Saybolt Viscosity Thermometer 
17F (66° to 80° F). The pressure gauge can be read to within + 0. 01 
psia and is accurate to within + per cent full scale. The thermome- 
ter can be read to within + 0.05° F. 

The flow measurement system consists of two subsystems 
of different flow measuring devices (Figure 2). The first is a paral- 
lel array of four rotameters, Een Models R-2-15-AA with stain- 
less steel float, R-2-25-D, R-2-25-A, and R-2-25-B with sapphire 
floats. These four rotameters are capable of measuring helium 
flow over the ranges 0.4 to 4.3, 1.0 to 10.9, 4.0 to 36, and 10 to 
108 cc/sec at 760 mm Hg and 70° Ee. The rotameters were calibrated 
by the author after they were installed in the facility by the water 


displacement method at atmospheric conditions. ach calibration 
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curve is based on at least 20 data points that are reproducible to 
within + 1 per cent from 30 to 100 per cent of full-scale reading. 
The overlapping range of adjacent rotameters permits each meter to 


be used in its upper, more accurate, range (Figure 7). 


The governing equation for a rotameter is: r 
2gV(p-p. ) 
e ee 
w 


where 
Q = volumetric rate of flow, 
V; = volume of float, 


g = gravitational constant, ; 


Pe = float density, 


ea fluid density 

fy = area of float, 

C = discharge coefficient, 

AL = area of annular orifice, = [(Dtby)*-d"] ; 

D = effective diameter of tube depending on position of 


float, 
b = change in tube diameter per unit change in height, 
d = maximum diameter of float, and 


height of float above zero position. 


wy 
I 


The value of C is a weak function of Reynolds number, but through 
careful design, the effect is largely eliminated for certain ranges of 
Reynolds number. Thus, when the roctameter is calibrated and used 


for a single fluid at approximately the same conditions of temperature 


' f 
® 
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and pressure, C may be taken as a constant. 
For a gas Pres is negligible compared with Pe» 80 equation 
(1) becomes 


ol 78 she 
Q — AL Aw e : (2) 
WwW 


Therefore, for any given position of ye float Q= Kp,-? where k is 
a ceneenen 

Once the meter is calibrated at known conditions of tem- 
perature and pressure for a given gas, the reading can be interpreted 
at other (not greatly different) conditions for the same gas by know- 


ing the gas densities or equation of state: 


Pwe 2 x r\ | 
a= a( e = Q (¢ z) (3) 


where the subscript c stands for the calibration condition. The 





mass flow rate (m) therefore can easily be determined: 


i 


m= Qp, = Q (oP! (2) 
The second flow measurement subsystem is a Brooks 
Vol-U-Meter. This is a precision-bore, vertical borosilicate glass 
tube containing a polyvinylchloride floating piston having a mercury 
O-ring seal. The usable length of the tube is 30 in. and the internal 
volume is indicated every 1 cc from Oto 25 cc. The precision of 
the instrument is 0. 2 per cent of indicated velene: By measuring 
the time required for the piston to rise through a given displace- 


ment and knowing the output pressure and temperature, the mass 


flow rate can be calculated. Volume flow rates from about 0. 003 
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to 2.5 cc/sec can be measured, the lower end being limited only by 
the patience of the operator and the overall stability of the system, | 
and the upper end by the reaction time of the operator anc the diffi- 
culty in gauging when the piston passes a given volume mark. 

It was found that the test gas stays at room temperature 
throughout the facility after it has pas sed through the water bath at 
room temperature. Thus, the temperature at the input to the flow 
measuring system is taken as the temperature in the Vol-U-Meter. 
A pressure tap is installed in the output line of the Vol-U-Meter. 
This pressure is indicated on a 0 to 25- and 25 to 50-psig dual- 
range dial gauge, Wallace and Tiernan Model FA2331i11, andis 
taken as the pressure in the Vol-U-Meter. It was found that the 
pressure drop across the mercury-sealed piston was about 0. 08 psi, 
or less than 1 per cent of the input pressure. 

During runs in which the Vol-U-Meter is to be used, the 
Vol-U-Meter is bypassed until a steady-state condition is attained. 
Once steady-state is established, the Vol-U-Meter bypass is closed 
and the test gas is drawn from the Vol-U-Meter until the flow rate 
has been determined. Once this has been accomplished, the bypass 
is again opened. The steady-state conditions are checked before, 
during, and after the volumetric measurement in order to be certain 
that the stagnation chamber remains at the same steady-state condi- 
tion. 

irom the flow measurement system, the test gas is fed to 
the stagnation chamber through 4 variable leak, Granvilie-Phillips 


series 203. This variable leak is used to controi the test-gas flow 
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rate, and thus the pressure in the stagnation chambex. The con- 
ductance of the variable leak is continuously variable from approxi-~ 
mately 100 to fom cc/sec with atmospheric pressure on the inlet 
side and vacuum on the outlet. Fully-closed conductance is stated 
by the manufacturer to be less than for cc/sec. Twenty-seven 
full turns of the driver handle provides precise control over the 
entire conductance range. A counter, graduated in units correspond- 
ing to 1/10 turn of the handle, is helpful in setting approximate con- 
ductances. At the higher flow rates and stagnation chamber pres- 
sures, a given setting gives reproducible performance. At the lower 
settings, hysteresis losses found in all devices with metal-to-metal 
seals cause slight variations in performance at any given setting. 
Figure 8 shows approximate counter setting versus stagnation cham- 


ber pressure Over the range of 0. U1 to Z00 mm He. 


Test Profile Plate 

mie test profiles are machined in steel Pietes which are 
fitted with O-ring seals into a cut-out in the dividing bulkhead of the 
vacuum vessel. The stagnation chamber side of the plates is flush 
‘ with the surface of the dividing bulkhead so that no undue entrance 
effects are present. The minimum diameter of the test profile is 
approximately 1mm. Thus, the ratio of tank diameter to profile 
diameter is about 400 tol. For all practical purposes, the stag- 
nation and exhaust chambers are of infinite size in comparison with 
the test profile. Therefore, the effects of the tank waiis on the flow 


field through the test profile are negligible. 


SO 
The first profile tested was a sharp-edged orifice witha 
diameter of 0.039 in. and an edge thickness of 0.0035 in. (Figure 9). 
Subsequent profiles will be short tubes and conical profiles of vari- 


ous entrance and exit half angles. 


Pump System 


iiie o-Inweourlet or thevexhaust Chamber 16 connected to 2 
welded aluminum tee. The tee has 4- and 2 1/4-in. flanged outlets. 
Both the outlets are fitted with handwheel-operated gate valves, 
Temescal Series 5000. From the 4-in. valve the flow drops 12 in. 
in a flanged Zee connector which serves as a baffle to prevent back- 
flow of oil from the Stokes pump. The Stokes pump is vibration iso- 
lated from the Zee connector by means of a 6-in. stainless steel bel- 
lows. The exhaust of the Stokes pump is vented outside the building 
(Figure 1). 

he .Z 1/4-in. Temescal valve is connected by O-ring seals 
Hestne input flange of the CVO diffusion pump. The ee of this 
pump is fed to the Cenco rotary backing pump by means of 1l-in. 
diameter copper tubing that is isolated from the diffusion pump by 
means of a short piece of Tygon tubing (Figure !). 

This pumping arrangement permits rapid pump-down of ° 
the system to the micron range by use of the Stokes pump. The 
pump-down time to 10 uy is approximately one minute (Figure 10). 
The diffusion pump can be isolated from the system and warmed to 
operating conditions while the Stokes pump is exhausting the system. 


Lhe switch-over to the diffusion pump is accomplished without any 
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increase in system pressure by simply opening the 2 1/4-in. Tem- 
escal valve and closing the 4-in. valve. It has been <ne practice to 
keep the diffusion pump in continuous operation and simply isolating 
it from the system when the flow rate and exhaust chamber pressures 
are such as to dictate the use of the Stokes pump as the primary 


vacuum source. 


Vacuum Pressure Measurement 

The primary means of measuring vacuum pressure is a 
bank of two McLeod gauges, made by Greiner Glassblowing Labora- 
tory, a O- to 800-mm mercury pence Wallace and Tiernan 
Model FA 135, and a 0- to 50-mm dial gauge, Wallace and Tiernan 
Model FA 160. One McLeod gauge is a non-linear gauge that meas- 
ures the range 0 to 125, ona 230-mm scale and that has a capture 
volume of 330.6 cc. The second, linear gauge, has three 500-mm 
scales that measure the ranges 0 tol, 0to 10, and 0 to 100 mm 
mataead Capture volume Of 250. 5 cc. 

In the non-linear gauge, a known volume of sample gas at 
the unknown pressure is captured and then compressed in a .-mm 
capillary that is sealed at its upper end until a reference column of 
mercury in a similar capillary is even with the end of the sealed 
Gapiiiary. lhe difference in the heights of the mercucy CCl mans 
the two capillaries is then a measure of the unknown pressure. As- 
suming isothermal compression, the following holds: 

PV = ov (5) 


where 





willie 


P = unknown pressure, 
V = capture volume, 
p = pressure in the compressed volume, 
“vy = compressed volume; and 
eee é (3) 
where 
A = capillary cross-sectional area, and 
h = difference in height of mercury in the two capillaries. 
With 
p = Pith (all pressures in mm Hg}, | (7) 
5 
eee eee 
V V 


from equations (5) - (7). 


In the micron range, P <<h , so that equation (38} becomes 


2 | 
_ Ah Z Z 
where k = known constant. 


In the linear gauge, a known volume is captured and then 
compressed to one of three known volumes. The height of a column 
of mercury above the compressed reference volume is then a meas- 
ure of the original unknown pressure. Assuming the isothermal re- 
lation again, PV = pv ; however, here v is the known reference vol- 
ume; then 


P = (Pth)< (10) 


er a 





ah 2 = 


eee (12) 
where k is a known constant. 

CVC GTC 100 thermistor gauges are used to monitor the 
stagnation and exhaust chamber pressures. These are used to indi- 
cate pressure variations in the tank and to check on the attainment of 
steady conditions, but are not used as a primary means of pressure 


measurement. 
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Lie ae ORY 

The flow of a gas through orifices, short tubes, and noz- 
zles at high pressure ratios in the slip anc transition regimes in- 
volves spanning the gap between gas cCynamics and gas kinetics. At 
both ends of this spectrum, the governing equations and boundary 
conditions are well defined, although solutions to the flow problem 
may not be readily obtainable. For Bite flow, the governing equations 
are the Navier-Stokes equations, but the boundary conditions are not 
known. The conditions at the wall of an object are no longer the zero 
svelocity Conditions of continuum flow. in the transition regime, the 
governing equations would seem to be those of gas kinetics, but ac- 
count must be taken of intermolecular as well as molecule-suriace 


collisions. 


Continuum Flow 

Far into the continuum regime where boundary layers oc- 
cupy only a small portion of the flow field, the flow problem can be 
analyzed to a good degree of approximation by assuming an‘inviscid 
Mewmwiield. in thisrepion, there exists a definite critical pressure 
ratio (>. /P .) beyond which the downstream conditions cease to in- 
fluence the upstream conditions. The maximum mass ilow is reached 
at this critical pressure ratio for givenwupstfeam Concitions« Ole 
smooth nozzle, this condition is reached wren the Mach number is 
unity at the throat, since oe maximum mass flow censity occurs at 
local sonic velocity. This’ maxiz. yressure ratio and correspond- 


. ~ seit - cc 
3g ifmaSS Ti0W are. 





ede 





Y 
e Y- 
Cc Yt1 
p= CP) (13) 
1 ytl 
P 2y-l 
2 
m = pAV = - + pA A 
\(T \/ you 
are 
i 5 ee 
Z 2 -1 
= p A(yRT.) a " ° (14) 


Although this choked condition gives the maximum mass flow, it does 
not determine the nature of the flow downstream of the throat which 
is essential to optimization of a nozzle. For this analysis, the 
methods of one-dimensional isentropic flow or characteristics can be 
used. In these analyses, the boundary layers on the walls can usu- 
ally be neglected to a good approximation. When they are taken into 
account, it is usually as a small correction to the nozzle area based / 
upon the boundary-layer displacement thickness. 

As the slip regime is approached, the boundary layers 
play a more important role in determining the over-all flow. Even- 
tually they fill the entire nozzle, and the gas flow is more like fully- 
developed laminar flow in a long tube than conventional nozzle flow. 
Marble has done a preliminary analysis which 

"attempts to treat the low Reynolds number nozzle prob- 
lem by modifying the formulation to include the viscous 
stresses associated with wall shear, but neglecting the 
normal viscous stress associated with acceleration 
along the flow direction. It is assumed that the flow is 
dissipative, but locally adiabatic. Under these restric- 
tions the results should be adequate to indicate the order 


of magnitude of the losses and of their effects on nozzle 
geometry, '! @ | 
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In this analysis, the velocity dependence on axial distance 
along the nozzle is chosen. This, teretier with the gas theered = 
namic properties and an estimate of the cross-sectional velocity 
profile, yields the information required to estimate the nozzle ge- 
ometry to achieve a given nozzle performance. Marble concludes 
that viscous dissipation plays a strong role for small Reynolds num- 
bers, and that for a given mass flow, as the length of the nozzle is 
increased, the area ratio must increase to achieve the same dis- 
charge velocity. For very low Reynolds numbers, the nozzle di- 
vergence angle must be increased as the length is increased. 

Although the choked condition for a nozzle in continuum 
flow occurs at modest pressure ratios (P_/P. = 1.89 for air and 
2.06 for He), much higher ratios are required for the flow through 
an orifice. This is because the sonic line for an orifice is S-shaped, 
while that of the nozzle is planar. Since the flow in the plane of the 
orifice is not sonic, .it follows that the mass flow is smaller than 
for a smooth nozzle of the same Reon area. 

The two-dimensional analogy to orifice flow, flow tirowee 
a slit, has been analyzed by the method of eam eetemietics the 
flow becomes sonic at the lip of the slit and flows locally around the 
lip like a Prandtl-Meyer expansion. The Mreamiine maps ona 
characteristic in the hodographic plane. When the characteristic 
connects the sonic point on the axis with the corner, the flow down- 
stream has no effect on the upstream conditions ana the flow is 


choked. This characteristic is an epicycloid from the sonic circle on 
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the u-axis and meets the corner characteristic at 9 = 45°. ‘he 
critical pressure ratio is:the one Perro eeondine to this 45° turn from 
sonic in Prandtl- Meyer flow. Liepmann gives the critical pressure 
ratios as 58.3, 25.6, and 21.5 for y's of 5/3, 7/5, and 4/3, re- 
peeciively. 
| Obtaining the mass flow is more difficult, for the full hodo- 
graph equation must be solved. Frankl has computed fora two- 
dimensional slit that the ratio of the mass flow to the mass flow ofa 
smooth nozzle of equal area is 0.85 for a gas an y=1.4. > 

The axisymmetric problem of the orifice is even more dif- 
ficult since the hodograph equations are not linear. Liepmann gives 
an argument that shows .that the required Prandtl- Meyer turning 
angle and pressure ratio are larger for an orifice while the mass 
flow is smaller than fora slit. The differences between the two 
should be small. Thus, it is expected that the mass flow through an 


orifice at high Reynolds number will be: 





1 yi 
3S 2 Y¥-l 
° oa Cc 2 
m = @ VT Vi fT A , (15) 
Cc 


where @ is on the order of 0. 85. ¢ 


As the Reynolds number is decreased, boundary layers 
build up near the lip and into the orifice.: These can have two effects. 
First, they reduce the effective area of ake orifice, and second, they 
can round out the orifice entrance and tend to increase the mass flow. 
Liepmann states thata limiting formula at high Reynoids number 


should have the form 
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el 
m = Mp ero! - CRe #) (16) 


‘where C is aconstant. He states that experimentally C: is of the 
order of 3. 

The third profile of interest, the short tube, would be ex- 
pected to act similarly to an orifice in the continuum regime with the 


Same entrance and boundary layer effects. 


Reece olecule Flow 

In the free-molécule regime, where the mean free path is 
much larger than the characteristic dimension of the profile in 
question, the flow through the profile can be determined by the meth- 
ods of kinetic theory. It is a matter of statistically keeping track of 
molecules and their surface interactions. 

Following Peotone one can compute the number of mole- 
_ cules passing through an orifice. Consider an element of area on the 
wall dA. Choose a system of spherical coordinates with origin at 
the center of dA and axis normal to dA. The mean free path of the 


gas molecules is 


KT 7 
=F (17) 


sy) 
V2 Pond” \/2 ae 
The mean molecular velocity is 
Cs. \/ (8/2) SS (18) 


The collision frequency per molecule is 


F =. V2 1a“cn (19) : 


d 








ag 

where n= molecule number density and a = collision diameter. The 
number of molecule collisions ina ee element dv eared ater, 
6, and 9¢ is equal to the average number of molecules in dv times 
the probability per unit time that ‘1 molecule collides, 

ndvF = (c/A)ndv . (20) 
Molecules leave dv in random directions. The fraction that leaves 
dv in the direction of dA is dA|cos6| [Ar“ - The fraction of mole- 
cules leaving dv toward dA which arrive at dA without undergoing 


further collisions is 


cndv dA| cos @ | e7tlh 


r fern 
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where e. is the probability that a molecule travels a distance r, 
from its last collision, without a further collision. The total number 


of molecules striking a unit area of wall per unit time is 


N == [fp teeseh om#! ay 7 (21) 
r 


4md 
ata WT 2 00 
= re | ae | d@ sin8 cos@ fa e7rlh 
0 0 0 
= = ; (22) 


Now, consider a small hole of diameter d much smaller 
than the mean free path. A molecule in the vicinity of the hole and 
traveling toward it will pass through the hole without suffering any 
collision in the neighborhood. The molecules pass through individu- | 


/ 


ally and independent of the motion of other molecules. 
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Equation (22) gives the number flow of molecules Eneouena 
a sharp-edged orifice when the rinieen number is sufficiently high. 
It can be applied to the gas on both sides of the orifice to give the net 
flow, since the two flows are ‘nd Seeedent of each other if the mean 
free paths on each side of the orifice are much larger than the diam-- 
eter. This type of flow is called effusive flow or effusion. If the gas 
on each side of the orifice is at the same tower eure, the net mass 


flow rate is therefore 


a A(P.-P.) (23) 


\/ 2nRT 


where 


= area of the orifice, 


upstream pressure, | : 


A 
A 

Cc 
ir 


e downstream pressure. 


When the orifice is not sharp-edged but has a finite length, 
some of the molecules collide with the side walls of the orifice and 
then are re-radiated in a random fashion. This causes some of the 
molecules that enter the orifice, or short tube, to be re-emitted to 
their original side of the orifice. Clausing showed that the come 


expression for circular cross-sections is 


ee) (24) 


\/ ome wu 


where x is a function of the length-to-diameter ratio of the orifice 

or tube and is called the Clausing factor. As the tube becomes long, 
7 : ae 

% approaches 4d/32 as anupper limit. Present gives a derivation 


for the mass flow for long tubes and shows that x= 4d/34 from geo- 
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metric and gas kinetic considerations . 
For long tubes of varying cross section A and perimeter H, 


the fundamental relation deduced by Knudsen is (cs. refs. 7, 8) 


ew Cc Oar (25) 
ie Oe ti‘ 
= dx 
A 
0 


For a conical shape of minimum diameter d, half angle @, and axi- 


allength 4, this relation gives 


l 
1 2 
ed Sr @ 1 ,d (d/2+Ltang) 
m = 7 (P.-P Mar) (aga) “dFttane) (26) 


The flow for other shapes in the free-molecule regime can 
likewise be calculated from geometric and gas kinetic considerations 


(cs. refs. 9,10). 


Slip and Transition Flow 


For flow in the slip and transition regimes, there exist a few 
formulae to calculate the flow through orifices and tubes. These are 
in the form of empirical or semi-empirical corrections to the con- 
tinuum equations. As an example, consider Poiseuille flow ofa gas 
through a long, straight pipe under the restrictions that (a) the gas is 
incompressible; (b) the flow is fully developed, i.e., the cross- 
sectional velocity profile is constant throughout the length; (c) there 
is no turbulent motion; and (d) the flow velocity at the tube wall is 
zero. The gas effectively flows in concentric cylinders which slide 
over each other. Since the flow is steady, the viscous shearing drag 


on the curved surface is balanced by the pressure. Therefore, 


~ 





ee 


2nrdx UL ca = -~wr° dP ° —— ca) 


With the velocity zero at the walls, integration of equation (27) gives 


2 d*-x® dP | : (28) 
“= Tou a | 
The total molecular flow (N) through the tube is | | 
d/2 4 
: eee 
N = n f u2nrar = 728 uKT es Aa oe 
0 ° aT 
and 
rc PS-PS 3 
(N = T38uKTt Z ) Zu 


since N is independent of.x by conservation of molecules. The mass 


flow rate is therefore 


of 
el 


ad 


— AP 
in Tees 


where P = (P_+P_)/2 and AP = Pee ; 

Although this Poiseuille equation is strictly limited by re- 
strictions given above, it has been modified for use in slip flow. The 
argument for this is given below. However, since the correction is 
based on experimental data, it should be pointed out that the modified 
equation must include a correction for compressibility effects, which 
are present, even though the argument is based solely on slip. | 

When the mean free path is smaller than the characteristic di- 
mension, but not negligibly small, an additive term to represent the 
slip at the surface is used. It is assumed that the gas at the wall has 
a finite velocity uo. The number of molecules striking a unit area of 


wall surface per unit time is 
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a? \ ap = 
N= nluot (rer) ae. Poe 


An estimate is made for Uo from mean free path considerations. It 


is found to be 


2k (au Zh da dp 
Tome) 95 ar) ~ “3 44 dx we 


where r = cylindrical radial coordinate from the axis of the tube, and 
x = cylindrical coordinate along the tube. The modified Poiseuille 


flow equation is (cf. ref. 6) 


. . wd = AP, me AP oe 
m ~ T28uRT 2 cu. 


The additive term has the same form as the free-molecule mass flow 
through a tube, but is smaller yr a factor ae Equation (33) is 
valid over a limited range and fails when the mean free path and 
characteristic dimension are approximately equal. 

‘, Dushman gives slip flow formulae for tubes attributed to’ 


a 
three sources. From Kennard he gives 








om 
SIP. nada —~ AP ae 
amen RT. (1+ > dhe Se 
where ¢€ is the coefficient of slip given by 
1 
SS ee 
_ a-f a 6 
re = 9 128) 


and f is the transfer ratio of momentum and represents the fraction 
of molecular collisions with the walls which result in diffuse scat- 
tering of the molecules. The value of f varies between about 0. 8 


and 1.0. 





IDB ice 
Knudsen found that the flow through a long tube can be 


written 


m =m, + 2m, | | (36) 


where ™m. is the continuum flow rate and my is the free molecule 
value. Z is a function of the mean pressure in the tube, the tem- 


perature, the tube radius, and the viscosity of the gas. 


1/ 
143 (gh) P 


Zs eens eeeta cae (37) 


1\l/25 
L+ 1.2355 (pe)? 


Sreekanth has recently developed a modified Poiseuille 
equation from experimental data on the flow through short tubes (/d 
from 0. 005 to 0. 995) at pressure ratios from approximately 1:1 to 
30:1 and at Knudsen numbers from 0.1 to 2 (the free-molecule end of 


the transition reemnen ce This modified equation is 


0.519 a? i 


AP) UT + a7e) (38) 


li = 
m = iene + = 


which is almost identical to the Poiseuille slip flow relation, equation 
(33), divided by the factor (1 + d/4). 

The above equations illustrate some of the empirical rela- 
tions that cover flows in parts of the slip and transition regimes. No 
known empirical relation or theory covers the entire range from 


free-~molecule to continuum flow. ' 
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IV. RESULTS 

Experiments were conducted to determine the mass flow 
through an Britice with an 4/d of 0.0895 (Figure 9) at Knudsen 
-numbers of from 8. 2 to 0.00068. The pressure ratios across the 
orifice were at ans 100:1. The experiments were conducted to 
verify that the low-density gas dynamic facility was functioning prop- 
erly and that the range from continuum to free-molecule flow could 
be covered. 

The initial set of data, which took one month to collect, 
mavenied that there was a slight leak into the stagnation chamber. 
The leak rate was approximately one micron-liter per second. This 
caused the measured flow rates to be too low near and in the free- 
molecule regime. The measured flow rates dropped to about 85 to 
90 per cent of the expected value. This small leak had virtually no 
effect on the data collected in the continuum regime. 

The initial test also revealed some difficulties with the 
flow measurement system and over-all system stabilization. When 
the Vol-U-Meter was used for flow measurement, the 1] per cent 
pressure change in the input to the variable leak when the flow was 
passed through the Vol-U-meter caused an upset in the steady-state 
previously ered This was remedied by throttling the gas in the 
feed line in parallel with the Vol-U-Meter to match the expected 
output pressure of the Vol-U-Meter. This meant that no change in 
feed pressure occurred when the flow was Gixected through the Vol- 


U-Meter. 
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At high Knudsen numbers it was found that it takes as 
much as five hours for the system to reach steady-state. This is 
due to the relatively large volumes and small flow rates. After the 
system stabilizes, it takes about one hour to collect one data point. 
The system is sensitive to temperature changes, the maximum al- 
lowable rate being about 2°F /hr. 

As the flow rate increases with decreasing Knudsen num- 
_ ber, the time constant of the system decreases and steady-state is 
attained fairly rapidly -- in as little as fifteen minutes. 

After the initial difficulties were remedied, anothe? series 
of data collecting runs were made. A typical test run is outlined in 
Appendix 1. Itis these data that will be discussed in the following 


sections. 


Data Reduction 

Two methods were used to reduce the data so that the 
results could be compared with the results of neo and 
Snaeen nee. 


Liepmann defines a characteristic velocity, W ; Mach 


number, Ma ; and Reynolds number, Re ; by 


ae + 
P-P i 
Ma= (EF tay 
P./e,)? c 


Re = —< = - (2) (41) 


26. 
where we is the kinematic viscosity of the gas evaluated at the up- 
stream conditions. He writes the mass flow rate in terms of the 


orifice area, A, and characteristic velocity, W, in the form 


ia Tp WA (42) 


where I is a dimensionless function which depends on: 


"Ma, Re, the ratio of specific heats y, and in 

the transition region possibly on slip, accomoda- 
tion coefficients, and possibly a Reynolds number 
based on the bulk viscosity coefficient. ''@ 


The parameter I is well defined in the two limits of con- 
tinuum and free-molecule flow. The free-molecule limit from — 
Knudsen's work is: 

Ma 
Te = aw . (43) 
In the continuum limit with Ma =1.0, 
yt} 4 


y-l 
Tr = ay) vo (44) 





where a(y) is theoretically well defined,varying from aq=1 fora 
. smooth nozzle to approximately a = 0.85 for an orifice. 
When the Clausing factor (x) is included to modify the | 


free-molecule limit, the relation for [ becomes 





ie = K ° | | (45) 


Liepmann presents his results in the form T/T 5. versus 1/Re. The 


,continuum limit for the former function is 


>= 


] L ve 


T/T, = - (2my)* (Ary — . (46). 


i 





apo. 
Sreekanth presents his results in the form m/m, 
versus Knudsen number defined as the ratio of mean free path evalu- 


ated at stagnation chamber conditions to diameter of the orifice or 


short tube. The expression used to calculate mean free path ae 


(84 )( 8 a a a 


As given by equation (24), the free-molecule mass flow limit is 


tel 


; it 
iie—— ae je 


This means that the continuum limit for the dependent variable is 
ous 
m a i, yo 
ow 2 Z ] 
= Sate a] i-Pye) = | ee) 


m 





In the limit with P. >>> P ,; this reduces to Liepmann's T/T 5. : 
It should be noted that the independent variables of Liep- 
mann and Sreekanth are closely related, and in fact differ by a small 


constant factor in the limit of ne >>> P. - The ratio Kn to mis is: 





lou. RT $ PP. 4 
Kn il d 
1/Re = (=p* so ak =) eee Po 4 “P uae oe i — 5. ¢-s4). (49) 


This reduces to 1. 275 when P . >>> Pp. 

Thus, the two systems of data reduction are similar, but 
not exactly the same. Both approach the same theoretical limits for 
ponent andere oleae flow. The present data were reduced 


by both of these schemes.. 


Data Presentation and Analysis 


Figure li shows the present data superimposed on Liep- 





iD as 
mann's. The two sets conform above 1/Re = Nou- and both seem to 
‘be slowly approaching the same free-molecule limit of 1.0. At the 
lower 1/Re end, the two sets seem to Be asymptotically approaching 
different raeentes This is as would be expected. The present data 
were for an orifice with 1/d = 0. 0895, while Liepmann's orifice had 
4/d = 0.0254. The Clausing factors for the two are, respectively, 
0.918 and 0.975. ‘ Assuming that a = 0.85, the two theoretical 
values for D/Ty. at the continuum limit are 1.68 and 1. 56, respec- 
tively. The last few values of T/T. obtained from the present data 
were equal to the theoretical limit of 1.68. Thus, both sets of data 
appear to be compatible at both limits and throughout the slip and 
transition regimes. They both show the same prior transition 
from continuum to free-molecule flow. 

Between 1/Re = 107¢ and 107+ , Liepmann's data exhibit 
a maximum which he attributes to the boundary layer smoothing the 
entrance condition to the orifice and thereby increasing the mass 
flow. The present data do not exhibit this maximum, but do exhibit 
a local maximum or a change in curvature which can be explained 
qualitatively by the same boundary layer argument. 

Figure 12 shows the present data together with Sreekanth's 
data for an orifice. The similarity between this presentation and 
Figure 1l is evident. Again, both sets of data are compatible. 
Sreekanth's data do not extend far enough into the low Knudsen num- 
ber range to get a good comparison throughout the slip and transition 

£ 


ranges, but a comparison at the free-molecwie end shows good agree- 


ment. 
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Sreekanth's empirical formula, equation (38), is also 
shown in Figure 12. Itis satisfactory at the high Knudsen soma: 
end of the transition regime. It overestimates the flow rate belowa 
Knudsen number of approximately 0.2. 

Over all, the present data show the transition from continu- 
um to free-molecule flow. The majority of the transition takes place 
within two orders of magnitude (0.01 < Kn<1, or i100<Re<1). 
Boundary layer and slip effects seem to be evident in the range 
Omw00l1 < Kn <0.01, or 1000 <Re < 100. The slow approach to the 


free-molecule limit extends beyond Kn = 10 or Re=0.1. 
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V. CONCLUSIONS 

| The continuous flow, low-density gas dynamic facility de- 
signed to study gas flow through various profile apertures in the 
slip and transition regimes was successfully designed, constructed, 
and tested. The test was performed with an orifice as the test pro- 
file. The data, which exhibit surprisingly little scatter, agree with 
previous data by Liepmann and Sreekanth. The facility can be used 
to gather data over a larger range than any previously reported. 
The range extends from the threshold of free-molecuie flow to con- 
tinuum flow. 

The collected data show that the transition to continuum 
from free-molecule flow takes place in essentially two orders of 
magnitude in Knudsen number. The transition is smooth, with no 
discontinuities as would be associated with an abrupt change in the 
flow field. A comparison of the data with the most recent empirical 
feerulae eee that there is no universal equation to cover the entire 
spectrum. The measurements show that the flow is approaching the 
theoretical limits at both high and low Knudsen numbers. 

Some interesting questions were raised by this study. The 
region around Kn = 0.01 or Re = 100 exhibits a local maximum or 
change in curvature that Liepmann reasoned was caused by boundary 
layers smoothing the orifice entrance conditions. This area is 
worthy of a more detailed investigation. The transition to free- 
molecule flow does not appear to be compiete ata Knudsen Nut oe csor 


iQ. The final approach to the theoretical free-molecuie limit appears 





aon = 
to occur very slowly. The final attainment of this theoretical limit 
also warrants further study. | 
In conclusion, it can be said that the facility ene its 
design criteria. It is hoped that it will be used to investigate the 
areas of orifice flow mentioned above and flows through other pro- 


files such as short tubes and nozzles. 
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APPENDIX IL. Typical Test Procedure 


I. Initial Pumpdown 


A. Commencing with all valves ciosed: 


I, Open shutoff valves to vacuum pressure measurement 
lines (3 valves) 
Z.. Open tank bypass valve 


— 


3. Open 4'' Temescal gate valve 


B. Start cooling water flow to: 


1, Stokes pump -- approximately 2 gal/min 


2. CVC diffusion pump -- approximately 1/15 gal/min 


ew mtart all pumps: 


U 


Gi 


a 


rh 


1, Stokes pump 
a. Check solenoid oil valve with ferromagnetic ma- 
terial 
b. If not magnetized, shut down immediately 
Z. Cenco Hypervac 25 backing pump 
3. CVC diffusion pump 
Fill McLeod gauge cold-trap reservoir with liquid nitrogen 
Open McLeod gauge stopcocks 
Close variable leak to a counter setting of 010. 
Open shutoff valve downstream of variable lonk 
When exhaust chamber pressure as read by thermistor no. 
is 10 microns: 


1. Close 4'' Temescal gate valve 


i 





ie 


is. 


Shut down Stokes pump 


I. Pump down system to desired ultimate pressure. Note: if 


this is the first pumpdown after the system has been ex- 


j « 1 LoL ; 
posed to atmospheric conditions 


dea 
ae 
Ne 


4. 


oe 


Open Stokes pump gas ballast valves (2) 

Delete steps D and £ 

Operate on gas ballast for 15 - 20 minutes 

Close gas ballast before going to Step H# 

Pump down system for at least 24 hours before pro- 
ceeding (it may be desirable to pump down for a few 


days to minimize outgassing) 


Execute steps D and E before proceeding 


il, Test Gas Flow Setup 


.A. Close tank bypass valve 


B. Close one of the shutofi valves to the McLeod gauge vacuum 


» 


-line (which one depends on whether exhaust or stagnation 


chamber pressure is to be measured) 


C. Fill feed line cold-trap reservoir with liquid nitrogen 


D. Start constant-temperature water bath 


He 


ee 


i. 


ey 


oe 


eemcuntenrrmmnee — ile 


Start motor and heater 


set temperature regulator after flow rate is established 


i. Open all valves in the gas feed line except: 


Close individual rotameter shutoff vaive 
Close rotameter cutlet vaive 
Close feed line shutofi valve immeciately downstream 


of the feed line vent vaive 





rr OO EC lm eee 


I, 
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Open helium (gas supply) bottle 
l. Set high-pressure regulator to about 100 »si (must be 
below 200 psi) 
Ze Set low-pressure regulator bo about Oe 
After a few seconds of venting feed gas to the atmosphere: 
1. Close feed-line vent valve 
2. Open feed-line shutoff valve downstream cf the vent 
Open variable leak to the desired counter reading (graph of 
approximate counter reading versus stagnation pressure is 
_ helpful for this) 
1. Monitor exhaust chamber pressure on thermistor no. l 
2. If exhaust chamber pressure exceeds 10 microns, 
switch pumps 
a. Close 2" Temescal gate valve 
-b. Start Stokes pump and wait 30 seconds 
c. Open 4" Temescal gate vaive 


When satisfied with pump setup, monitor 


stagnation cham- 


ber pressure on thermistor no. 2, 0 - 50 mm dial gauge, 
or manometer. Note: Keep manometer shutoff valve 
closed below 50 mm Hg to minimize m 


tarnination oi the variable leak. 


LI. Flow Rate Measurement 


A. Low Flow Rate -- Vol-U-Metez 


1, Open input valve to rotametesr no. 


€. Open rotameter output-line sEutocil velve 





~43~ 


3. Close rotameter bypass valve 


4. Partially open rotameter no. 1 needle valve to give 
| ' 
about .j psi pressure drop between flow measure- 


ment system input pressure (0 - 15 psia dial gauge) 


and Vol-U-Meter output pressure (0 - 25, 25 - 50 psia 


dial gauge) 


5. Stabilize system for at least 2 hours 


6. Continually refill cold traps 

7. Close Vol-U-Meter byvass valve anc simultaneously 

: open rotameter bypass valve. Note: this procedure 
minimizes changes in the gas input pressure to the 

,. variable leak when the Vol-U- Meter is thrown into the 
SYoteim. 

8. Time the piston between any two volume marks (pref- 
erably 0 to 25 cc at the slower rates, and O and 2Zcc 
at the faster rates) 

9. Before the piston reaches the top of the tube: 

a. Stop the piston by cracking open the Vol-U-Meter 
bypass valve. 
b. Lower the piston gently to the bottom of the tube by 


carefully opening the Vol-U-Metex oypass valve a 


bit more 


B, Moderate to High Flow Rates -- Rotamecvers 


i. Close input needle vaive of desired rotameter 


& Open input valve to rcotameter 





SAS 
3. Open rotameter output-line shutofi vaive 
4. Close rotameter bypass valve 
5. Open rotameter input needle valve carefully to full 
open position. Note: this should be done carefully to 
prevent damage to the rotameter float. 
the rotameter can be reac continuously and will help 
indicate when the file is stabilized (tre >rimary meens 
of assuring stabilization is stagnation chamber pres- 
sure measurements). 
IV. Other Data Measurements 
A. Atmospheric conditions 
1, Pressure - barometer (apoly temperature/gravity cor- 
rection to indicated Hg column height) 


Z. Lemperature - thermometer on barometer case 


B. Gas inlet conditions 
1, Pressure - before flow measurement system 

a Low-pressure regulator dial gauge (0 - 30 in. H,0) 
- this pressure plus atmospheric should equal flow 
measurement system input pressure 

b. Flow measurement systera input cial gauge (0 - 15 
psia) 

ce Flow measurement system input manometer (0 - 
800 mm Hg gauge) 


3. Temperature - before flow measurement eiisce m4 
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C. Stagnation Chamber Conditions 
1, Pressure 
ae YThermistor no. 2 


b. McLeod gauges 


c. Dial gauge (0 - 50 mm Hg} 


d. Manometer (0 - 800 mm Hg) - the scale is already 


, Gotrectecween temperate oases 


SCalGereacding trom coumecred ba 


ravity - subtract 


POS ener caaino | 


Zz. Temperature - thermometer is in protective shroud 


protruding from uvstream end of tank 


D. Exhaust Chamber Conditions - Pressure 
1. Thermistor no. 1 
2. McLeod gauges 


V. McLeod Gauge Operation 


A. Activate pressurization system 


1, Close individual gauge shutofi valves 


2. Open vent line valve 
3. Open needle control valve 
4. Open nitrogen bottie valve 


5. Open regulator shutoff valve 


6. Set regulator pressure at approximately 25 psig 


7 Close needle control vaive 
B. Sample measurement 
i, Capture of gas sample (either gau 


-a.e Close vent line valve 


b. Open desired gauge shutofi valve 


aeulie 


c. . Open needle valve slowly to obtain rate of mercury 
rise such that sample capture takes 15 - 30 sec- 


,onds 


‘2. Non-linear gauge 


a After sample capture, increase rate of mercury 
rise but slow it drastically near the top of the 
large volume to avoid shocking the closed capil- 
lary tube 

b. When mercury is in the closed capillavy, increase 

them@rave ot rise again 

Cc. wlow the mercury rise as fhe ievel in the open capil- 

lary nears the top 

d. Stop the top of the mercury meniscus at the bottom 
of the black tape line by either closing the needle 
valve or closing the gauge shutofi valve 

e. Tap both capillaries a few times to overcome sur- 
face tension effects 

fr Keaa the dimference in mercury Column Neients 

3. ‘Linear gauge 

a. After sample capture, increase rate of mercury 
rise but slow it near the top of each large volume 
to avoid overshooting the compression volume cal- 
ibration mark 


b. Stop the mercury meniscus at tHe des: seqvecrived 


. - 4 q S - =<. = @ dey 
volume marie@oy cle@sine cities tac 
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valve ar tho needle valve 

(1) Use Ist scribe mark omc leit tube for pres- 
sures between 125 end 10 mm 

{2) Use “Aa scribe mark and micdie tube for 
pressures between 10 and 1 mm 


_ (3) Use 3rd scribe mark and right tube for pres- 


sures between.! and.! mm (pressures below 


r 


- 15 mm can also be read on the non-linear 


() 


McLeod) 
G. Read the difference in mercury column Heres 
4, -To lower the mercury level in either gauge: 
a. Close needle valve 
b. Open vent line vaive 
c. Open gauge shutoff valve 
IN@te: 12 a slug of Mercury bens up i the closea capil- 
lary of the non-linear gauge, carefully neat the tube 
with a fuel-rich (yellow) natural gas /oxygen flame until 
the slug vaporizes. 
5. Conversion from scale reading to pressure 
a. Apply temperature and gra -ty correction to all 


mercury Column heights beiore using scalc saccore 


G.. Scale factors 


: ; -Z Z 
(1) 0 - 100 micren cauge® @, = . 2a eee 
- ~LeeS aw. 
(2) O-i ram gauge: P, , = Ont 
Racasans PPOs Bt 
Ca’ _—, - ~ Sep 2 -Z. 
(3) 0 = 10 mm pauce: FP, , =. S92 ek 
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(4) O-100 mr gauge: P = 2.494x10 “h 
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Close valve downstream of variable leak 


Close feed-line shutoff valve GCownstream of feed line 


vent 


Check Vol-U-Meter bypass and rotameter bypass 


e 
4 


valves open 
Turn off helium (test gas) bottle 


Open feed-line vent valve 


B. CVC Diffusion purp 


Close 2" Temescal gate valve 

Turn off diffusion-pump heater switch 

Let water and Hypervac 25 run until neater feels cool 
to the touch - 

Turn off water 


Lurn off Hypervac 25 


C. McLeod pressure system 


Close both McLeod gauge shutoff valves 
Oven vent line valve 
ilurn off nitrogen bottle 


Open needle control valve 


Vil. Periodic Checks 


y 2 


A. Check wooden shims under the test-texk support stand twice 


& 4 


I q —— % t x bY ~~ % . he — a a ee ; ms rs. yf ~ et 
. weekly and krock them back into positions markec on the 
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floor if vibration has moved them out. 


Check the oil level in the Hypervac 25 baciting pump at 


least every other day. Add Cenco HyVae oi. while pump is 
running. | 

Oil level in the Stokes 412H pump. shoulc be halfway up the 
sight glass while the pump is running. 

Temperature of diffusion-pump cooling water taxen at the 
outlet should be 110 - 120 degrees Fahrenheit for maxi- 
mum pumping rate. Check daily. 

Stokes pump cooling water should be 70 - i100 cerrees 
Pahrenheit at the outlet. Check cee eaer period.of ex- 
tended operation. 


Drain the water out of the Stckes pump and refill with clean 


water twice weekly. 
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